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ABSTRACT 

Satellite soundings of the topside ionosphere reveal diffuse signals on the 
low frequency side of nf H when n^ < f T where n is an integer and f H and f T are 
the electron cyclotron and upper hybrid frequencies, respectively. They extend 
no further than the mid-way point between (n - l)f H and nf H . The present observ- 
ing conditions allow detection up to n = 10. The origin of the signals is attributed 
to the Harris instability of longitudinal plasma waves when /T > 2n where 
Tj- and T 1( refer to the temperature corresponding to the electron motion per- 
pendicular and parallel to the earth's magnetic field, respectively. The large 
anisotropies are caused by collisionless cyclotron damping of energy from the 
high power sounder pulse. This energy is absorbed in a region which extends 
radially from the antenna a distance which is small compared to the antenna 
length but large compared to the antenna sheath region. 
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ELECTROSTATIC VELOCITY-SPACE INSTABILITIES STIMULATED 


NEAR THE HARMONICS OF THE ELECTRON CYCLOTRON 
FREQUENCY IN THE IONOSPHERE 

INTRODUCTION 

In 1959 Harris showed that obliquely traveling plasma waves, consisting of 
longitudinal oscillations of electrons, can be unstable in the presence of an ani- 
sotropic electron velocity distribution 1 . These instabilities occur near the 
harmonics of the electron cyclotron frequency ^ . This work has been extended 
to include ion motions with the result that instabilities are predicted at harmonics 
of the ion cyclotron frequency 2, 3 . The theory has been related to plasma con- 
finement experiments that produce large anisotropies in the ion velocity distri- 
bution in an attempt to explain the enhanced radiation often observed at these 
frequencies 4 . It is often difficult to interpret laboratory observations in terms 
of a single instability mechanism; the intense microwave radiation observed at 
nf (n = 1 to 19) by Wharton et al 5 is a case in point. It has been suggested that 
many of the laboratory observations may be explained in terms of a double dis- 
tribution, i.e., where a hot anisotropic component of ions (or electrons) are 
superimposed on a cooler nearly isotropic component of the same species 6 . 

This interpretation has been used by Schwartz and Lichtenberg 7 to explain 
their observations related to velocity-space instabilities at nf j( with n = 1, 2, and 
3. The present paper will present observations in the ionospheric plasma which 
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provide evidence for the stimulation of Harris-type instabilities at nf H in a single 
electron distribution. 

The original work of Harris 1 assumed Tj/T,| = « where Tj_ is the tempera- 
ture corresponding to electron motion perpendicular to the earth's magnetic field 
vector fi and T|j is the temperature corresponding to electron motion parallel to 

When finite T_l/T|| values are considered, the instability region extends be- 
low the nf H value to a point midway between nf H and (n - l)f H for short wave- 
lengths 8,9 . From the work of Gruber et al. 10 and Oya 11 it can be seen that the 
growth rate of the instability increases as the wavelength decreases; this be- 
havior is illustrated schematically in Fig. 1 where the normalized dispersion 
curves are presented (the wave number being normalized by the inverse of the 
electron cyclotron radius R). Thus the Harris type instability is excited most 
readily when the frequency is near the mid point between the harmonics of nf H . 
Also illustrated in Fig. 1 is the change in shape of the. dispersion curve for a 
given nf H as the upper hybrid frequency f T crosses the nf H value. It is evi- 
dent from the figure that the instability can exist near nf H only when f T > nf H . 

OBSERVATIONS 

The observations are from the ISIS 1 satellite, ISIS 1 is the third of 4 
satellites of the series designated as International Satellites for Ionospheric 
Studies (the other three are Alouette 1, Alouette 2, and ISIS 2). Each satellite 
contains a sweep frequency rf sounder in order to obtain electron density vs. 
altitude profiles in the topside ionosphere, i.e., the region above the electron 
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density maximum near 300 km. The later satellites contain additional radio 
experiments, probes, particle detectors, and optical experiments and have evolved 
into extremely sophisticated ionospheric laboratories . 12 Of particular interest 
to the present study is the sounder experiment because, in addition to the electro- 
magnetic waves used to obtain the electron density profiles, it generates plasma 

waves of low group velocity at fundamental frequencies in the ambient plasma. 
These received waves have been described as ''resonances" because 

they appear as long duration signals (many milliseconds) which are re- 
ceived immediately after the transmission of a short (100 m sec) rf pulse. An 

example of a swept -frequency recording (called an ionogram) that clearly illus- 
trates these resonances is presented in Fig. 2 . An ionogram is composed of a 
large number of individual vertical scan lines, each at a separate frequency. 
Resonances are observed at the electron plasma frequency f N , the upper hybrid 
frequency f T , and at nf H (n = 1 , 2, 3 * • * ) where, using the conventional notation 
from the ionosphere literature, f N (kHz) «s 9 [N(cm -3 )] l/l , f H (kHz) * (1/36) B( 7 ), 
f 2 = f 2 + f 2 , and N and B are the ambient values for the electron density and 
magnetic field strength respectively (I 7 = 10 -5 Gauss). The resonances at f N 
and f T have been interpreted in terms of the reception of obliquely traveling 
plasma waves of low group velocity which return to the sounder because of the 
local gradient of N, as first proposed by McAfee 13 ; the theory has received con- 
siderable observational support 14 . Andrews and Fang 15 suggested that the same 
oblique echo mechanism may also be applied to the nf H resonances (n > 1 ) whereas 
the ray tracing calculations of Muldrew 16 indicate that the gradient of f N is of 
secondary importance compared with the gradient of f H for these resonances. 
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The approach of Muldrew may account for the long durations observed with the 
low order harmonics of f H but such propagation effects may not be as important 
for the higher order harmonics. 17 Also shown on Fig. 2 are resonances labeled 
f mj f Q3 , and f Q5 . These resonances are members of sequences of resonances 
that appear between the harmonics of f H ; the f Q series 18 is observed above f T 
and the f D series 19 is observed below f T . The f Q resonances have been inter- 
preted in terms of stimulated plasma waves with group velocities that are very 
nearly equal to the satellite velocity , 20 whereas the f D resonances involve 
plasma waves with group velocities approximately an order of magnitude larger 
than the satellite velocity (which is typically 7 m/msec). The f D resonances 
have been interpreted by Oya 11 in terms of the instability introduced by Harris 1 
together with a nonlinear interaction of the plasma waves associated with the f D 
and f q sequences of resonances and the 2f H resonance. 

In order to illustrate the received signals that are attributed to the Harris 
instability, two ISIS 1 satellite passes were selected which contained data cover- 
ing a wide f N range and a narrow f range. This situation is desirable because 
it allows the spectral changes of a nearly fixed frequency pattern of nf H resonances 
to be observed as f T sweeps through the pattern. The variation of the plasma 
parameter ^,/f H for the data of interest is shown in Fig. 3 together with the 
variation of satellite altitude and the change in dip latitude. (The dip latitude is 
obtained from the dipole field expression tan" 1 (1/2 tan I) when I is set equal to 
the inclination of the true magnetic field direction from the horizontal.) The 
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corresponding ionospheric data are presented in Figures 4 and 5. Fig. 4 is a 
composite presentation of a series of ionograms corresponding to the range of 
fpj/fjj given in Fig. 3a; similarly, Fig. 5 corresponds to Fig. 3b. Together, the 

observations illustrate the change in spectral composition of each nf H resonance 

between n = 2 and n = 10 as the plasma conditions change from f > nf to f < nf . 

T H T H 

The f T resonance (designated by T in the figures) occurs between the plasma 
frequency (designated by P) and the exit frequency of the electromagnetic x trace 
(designated by X); it is not labeled on the ionograms where it is not a prominent 
feature. The main features to note in Figures 4 and 5 are the following: 

(1) The nijj resonances are observed over a wider frequency range and with 
a longer time duration when n^ < f y than when nf H > f T . 

(2) In the former case, a diffuse signal is observed to extend from nf to 

H 

lower frequencies but not below the mid-point between (n - l)f and nf . 

(3) In general, this diffuse signal appendage is observed to be strongest on 
the nf H resonances nearest to, but less than, f T . 

(4) The time duration of the diffuse signal increases as the frequency 
approaches nJ^ . 

(5) The nf H signal is greatly enhanced when the exit frequency of the 
electromagnetic z trace at the satellite coincides with nf , 23 

H 

The first two observational points agree with the theoretical expectations 
for the Harris instability as outlined in the introduction. The third point may be 

related to the calculations of Oya 11 which indicate a slight displacement of the 

instability region to higher kR values for a given anisotropy as f N /f H increases. 

A higher Tj_/T|| ratio would then be required to produce the instability near nf H 
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(low kR value) when f T > > nf H than when f T > nf H . One must use caution, 
however, when interpreting the observations for the dependence of signal 
strength on fj^/f H since instrumental effects can be important. 21 The fourth 
point appears opposite to expectations since the instability growth rate is 
expected to be greater for the frequency range near the midpoint between the 
nf, , values rather than near nf H (see Fig. 1). The capability of generating the 
required anisotropy must also be considered and, in the next section, it will be 
shown that the anisotropy produced by the transmitter is much greater near nf H . 
Also, the finite antenna size may play a fundamental role in determining the 
spectral response since the wavelength of the oscillations increases as the 
frequency approaches nf H « The fifth point implies that more energy can be fed 
into the electrostatic mode when it coincides in frequency with the electromagnetic 
cutoff of the z trace. This observation has been noted earlier 17 . 

DISCUSSION 

The main question pertaining to the interpretation of the present observations 
in terms of the Harris instability concerns the generation of the required elec- 
tron temperature anisotropy. The anisotropy required for the instability near 
nf H is given by 8 


V T II 


> 2n. 


( 1 ) 
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Large anisotropies of this type can be created by the rf electric field of the 
sounder pulse since more energy is absorbed by the component of electron mo- 
tion perpendicular to B than by the component parallel to B due to the collision- 
less cyclotron damping process. The magnitude of this damping is illustrated 
in Fig. 6. The absorption is very severe and it increases as the propagation 
angle 9 decreases (cyclotron damping is not present for propagation perpendicular 
to B, i.e., when 6 = 90°). The figure corresponds to k^R fixed and 6 variable, 

6 decreases as the frequency approaches the value of 2f H . This behavior is 
illustrated very vividly by the calculations of Oya 11 which show that the dispersion 
curves are much closer to the nf H frequency value when 8 << 90° than when 
9 = 90° for the condition of interest, i.e., f T > nf H . Thus the cyclotron damp- 
ing of plasma waves is much greater near the harmonics of f H than it is in the 
frequency region mid way between the harmonics. This increased cyclotron 
damping near nf H could explain the increased signal strength observed near nf H 
(as pointed out in item 4 of the last section) , since the amount of transmitted 
energy absorbed by the medium increases dramatically as the sounder frequency 
approaches nf H . The result of this energy absorption will be to increase Tj^ in 

the region around the sounder antenna. 

In order to estimate whether this increase in T^ is significant with respect 
to the thermal electron energy, it is necessary to know the effective plasma 
volume for energy absorption and the energy of the transmitted wave . The 
magnitude of the absorption indicated in Fig. 6 (i.e., 10’ s of dB per cyclotron 
period) suggests that the effective volume of enhanced T^, due to the 
absorption of plasma waves propagating with low group velocity, is very small. 
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Observations 17 indicate that the signals responsible for the higher order 
nf resonances have near zero group velocity and that they do not extend beyond 
about two meters in the radial direction from the antenna element. These con- 
clusions were based on the nf^ resonances excited whennf H > f T , i.e., when the 
instability condition is not satisfied. The present observations indicate that the 
nf H resonances persist for a longer time duration when n^ < f T due to the Harris 
instability. If, however, it is assumed that the instability region extends over a 
region of space of the same order of magnitude in size as the region covered by 
the plasma waves when no instability is present, then an estimate of the energy 
density absorbed from the rf sounder pulse can be made. Next, it is necessary 
to make an estimate of the fraction of the transmitted energy that wall be available 
in the electrostatic mode (which will be the effective mode for energy absorption 
in the region near the antenna). It is apparent from the regular reception of 
ionospheric echoes over distances of the order of thousands of kilometers (e.g., 
see Figures 4 & 5) that a very significant amount of energy is radiated in the 
electromagnetic mode. It is also apparent from the intensity of the plasma 
resonance signals (again - see Figures 4 & 5) that a very significant amount of 
energy is radiated in the electrostatic mode. (Evidence for significant energy 
radiation in the form of electrostatic waves from an antenna immersed in the 
ionospheric plasma is also available from in situ impedance measurements - as 
first pointed out by Whale 25 ). For the sake of discussion, assume that 1/2 of the 
total energy contained in the sounder pulse is available in the electrostatic mode 
and that this energy is absorbed in a cylindrical volume extending a radial 
distance r from the antenna of tip-to-tip length L. Then the absorbed energy 
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density is (P/2) 2 L) where the ISIS 1 sounder power P and pulse length T 

are 400 w and 98 /xsec, respectively. When r = 2m and L = 73m (for the long 
dipole antenna), the absorbed energy density is 2 x 10 “ 4 ergs/cm 3 which is 
more than 3 orders of magnitude larger than the thermal energy density of 10 7 
ergs/cm 3 based on N = 3.6 x 10 5 cm -3 and T = 2000°K (as appropriate for the 
data of Figure 4d). Such a large value for the energy absorption indicates that 
the estimates for the fraction of transmitted power and for the size of the 
absorption region have considerable leeway while the instability anisotropy 
requirement of (1) remains satisfied. 

The above estimate of a few meters for the region of enhanced T^/T|| is 
comparable to the distance covered by the wave packet propagating in the elec- 
trostatic mode during a time interval equal to the sounder pulse duration 26 . 

This time interval is equal to about 60 times the electron cyclotron period and 
the results of Fig. 6 indicate that waves propagating at oblique angles with 
d« 90° would be completely damped out in a much shorter time. The calcula- 
tions leading to Fig. 6, however, are appropriate only for a particular value 
of k^R (namely k x R = 1 . 0) and for an isotropic plasma. When smaller 
kjR values are considered, then the dispersion curves corresponding to 9 < 90° 
(where the least cyclotron damping occurs) approach the frequency at nf H u . 
Also, in an anisotropic plasma with Tj_/T|| > 1 the cyclotron damping is reduced 
as can be seen by writing the condition for negligible cyclotron damping 


c«j + n« H 

/ \ 1/2 
f m ) 

k !l 

\ 2 kT ii/ 
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in the form: 


k„R < < 


v t iiV /2 i« 


where k is Boltzmann's constant, to = 2irf, cj h = 2irf H , and R = (fcT^/m) 1,4 co^ 1 
is the cyclotron radius. Using (2), and the dispersion equation solutions for 
Tj^/Tj| = 1 and X ± /Tj| = 10 as given in Figures 6 and 8b of Oya 11 , it can be shown 
that the range of propagation angles satisfying the condition of negligible cyclo- 
tron damping increases with increasing Tj/Ty . This result indicates that as 
energy is absorbed from the leading portion of the rf pulse , an anisotropy is 
produced which decreases the cyclotron damping for the later portions of the 
pulse which can then propagate further into the medium. This reduction of cy- 
clotron damping with an increase in T^/Ty is not very strong, however, because 
of the power of l/2 in (2). Thus the figure of a few meters for the extent of the 
region of large Tj_/Ty is a reasonable upper limit, especially near nf H where 
the damping is extreme for certain combinations of kj_R and 9. 

The above discussion suggests that the temperature anisotropies required 
to produce the Harris instability could be produced by the cyclotron damping of 
electrostatic cyclotron harmonic waves in the vicinity of the sounder antenna, 
i.e., within a radial distance of a few meters from the antenna element for the 
higher order harmonics of f H . It should be noted that this distance is well outside 
the antenna sheath region which is of the order of a few centimeters,, The diffuse 
signals associated with the nfn resonances (when nf H < f T ) in Figures 4 and 5 


are attributed to this instability. 
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Figure 1. Schematic normalized dispersion curves for the electron electrostatic cyclotron harmonic waves. 
The dashed curves correspond to waves that are not subject to the Harris instability whereas the solid curves 
are subject to this instability if the required electron temperature anisotropy is present. The increasing 
thickness of the solid curves indicates the increasing instability growth rate with increasing wave number (de- 
creasing wavelength). Note the change in the location of the dispersion curve associated with (n +■ 1) f H from 
the stable region above (n + 1) f H when f T < (n + 1) f H in (a) to the unstable region below (n + 1) fT when 
f T > (n + 1) f H in (b). 
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Figure 2. The high resolution portion of an Alouette 2 ionogram recorded at the Quito Telemetry Sta 
tion on 4 May 1967 (12:14:22 UT; 25. 2°S, 65. 8°W, 2612 km in altitude). The heavy vertical lines are 
frequency markers and the horizontal lines are range markers separated by 200 km in range or 4/3 
msec of delay time. The heavy vertical traces, which are due to electrostatic waves of low group 
velocity in the vicinity of the sounder, are identified at the top of the ionogram; the weaker non ver- 
tical traces are due to the ionospheric reflection of electromagnetic waves radiated by the sounder. 
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Figure 3. The variation of the ambient plasma parameter f N /f H , the ISIS 1 altitude, and the dip latitude as a 
function of universal time for (a) QUI pass 1776 on 7 July 1969 and (b) SNT pass 2033 on 30 July 1969. The 
letters a-f on the f N /f H curves correspond to the ionogram data presented in Figures 4 and 5. (The dashed 
portion of the f N /f H curve in (a) corresponds to a period of missing data due to interference, ) 




Figure 4. A series of ISIS 1 ionograms from QUI pass 1776 corresponding to 
the conditions given in Figure 3a. The letters 1 through 14 are used to identify 
the resonances at f H through 14 f H , P represents the plasma frequency reson- 
ance at f N , Z and X represent the cut-off frequency at the sateHiteforthee^ec- 
tromagnetic z and x traces, respectively. Certain instrumental effects 
must be considered when evaluating these ionograms. 
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.Figure 5. A series of ISIS 1 ionograms from SNT pass 2033 corresponding to 
the conditions given in Figure 3b. Dj and D 2 represent the diffuse resonances 
at f D1 and f D2 , respectively. Q 3 represents the electrostatic resonance at fq 3 
and T represents the upper hybrid resonances at f T . 
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Figure 6. The plasma wave absorption A due to cyclotron damping as a function 
of f/f H for the case where (f N /f H ) 2 = 5. 0 and k^R = 1.0 (the propagation angle 
decreases from 79° to 63° as f/f H increases from 1. 5 to 1. 86). (Adapted from 
the calculations of Tataronis and Crawford i4 .) 





